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Recent theoretical and experimental work has demonstrated the existence of one-sided, invariant barriers to 
the propagation of reaction-diffusion fronts in quasi- two-dimensional periodically-driven fluid flows. These 
barriers were called burning invariant manifolds (BIMs). We provide a detailed theoretical analysis of BIMs, 
providing criteria for their existence, a classification of their stability, a formalization of their barrier property, 
and mechanisms by which the barriers can be circumvented. This analysis assumes the sharp front limit and 
negligible feedback of the front on the fluid velocity. A low-dimensional dynamical systems analysis provides 
the core of our results. 



The passive advection of particles in a two- 
dimensional fluid flow is a well-studied prob- 
lem important both for its direct relevance to 
fluid mixing and for its broader connection to 
Hamiltonian dynamics. The latter finds applica- 
tions ranging from chaotic ionization in atomic 
physics^J to the motion of asteroids in the so- 
lar syste A key observation in passive chaotic 
advection is that particular curves — the invari- 
ant manifolds — form barriers to particle trans- 
port 9 12 . The question remains to what extent 
an active material within a fluid flow is subject to 
such barriers. We are particularly concerned here 
with active systems that generate fronts, e.g., 
chemical react ions^ffill or sound waves in moving 
fluids, plankton bloomd^E^ in ocean flows, opti- 
cal pulses in an evolving optical medium, phase 
transitions in liquid crystals, or even the spread- 
ing of disease within a mixing population. Here, 
we adopt a model for the propagation of points 
along the front in which each point is viewed as an 
independent active "swimming" rod — while the 
rod always swims perpendicular to its orienta- 
tion, the fluid current both translates and reori- 
ents the rod. The invariant manifolds that nat- 
urally arise for this dynamics are a modification 
of the traditional manifolds governing passive ad- 
vection. They prove to be fundamental, one-sided 
barriers to front propagation and strongly influ- 
ence the patterns of front evolution. 



I. INTRODUCTION 

Recently, the existence of robust barriers to the prop- 
agation of reaction fronts in advection-reaction-diffusion 
systems was experimentally demonstrated in both time- 
independent and time-periodic vortex-dominated flows^. 
The flows were magneto-hydrodynamically-generated, 
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quasi-two-dimensional, and vortex-dominated, on cen- 
timeter length scales. The reaction was a ferroin- 
catalyzed, Belouzov-Zhabotinsk}B^ reaction. See the 
accompanying article by Bargteil and SolomorP^ for 
additional experimental evidence in irregular, vortex- 
dominated flows. In Ref. 18, the barriers were experi- 
mentally shown to be both (i) one-sided, preventing front 
propagation in one direction, but not the reverse and (ii) 
invariant, staying either fixed in the lab frame or period- 
ically varying (for periodically driven flows). These bar- 
riers therefore play an important role in pattern forma- 
tion in such reactive systems; they also limit, and could 
potentially be used to control, the average reaction prop- 
agation speed through the fluid. 

The explanation for the barriers in Ref. 18 was based 
on a low-dimensional dynamical systems approach, which 
interpreted the experimental barriers as invariant man- 
ifolds attached to unstable periodic orbits. Hence, 
these barriers were named burning invariant manifolds 
(BIMs), where "burning"^ emphasizes their relevance to 
reaction-front propagation and distinguishes them from 
the well-established invariant manifolds that act as bar- 
riers in the advective transport of passive tracers 9 ^ The 
main objective of the present paper is to develop a ped- 
agogical and comprehensive theoretical justification for, 
and elaboration of, the picture introduced in Ref. [18] . In 
particular, we clarify the physical properties of BIMs as 
bounding curves and give explicit criteria (in the time- 
independent case) for the existence of BIMs in laboratory 
flows. 

A complete model of a general advection-reaction- 
diffusion system is typically built, at least initially, 
around a partial differential equation that couples the 
local chemical dynamics to the local fluid flow and vice 
versa. In general, the reaction dynamics can profoundly 
impact the fluid velocity (e.g. in combustion.) However, 
such feedback appears to be small for the expe rimen ts, 
such as Refs. 18 and 20, motivating this worl J 19 ' 22 ^ . 
Thus, in this paper, we take the fluid velocity field u(r, t) 
to be prescribed at the outset and to be unaffected by the 
subsequent chemical dynamics. Furthermore, we take the 
"geometric optics" limrP^ in which the reaction is suffi- 
ciently fast that a sharp front develops between those 
regions that are fully reacted and those that are com- 
pletely unreacted. These assumptions lead to a three- 




FIG. 1. (color online) a) Fronts converging (blue to green) 
upon BIMs (bold, red) in a time-independent vortex flow. 
Panels b-d show fronts initialized at different stimulation 
points (small open circles) burning through BIMs opposite 
their burning direction but bounded by BIMs in their burn- 
ing direction. Panel d includes a BIM attached to a burning 
fixed point near the top of the channel, which restricts the 
front in the upper right. 



dimensional ODE for the propagation of a point on the 
front, discussed below. 

For passive tracers, it has long been recognized 
that one-dimensional invariant manifolds bound advec- 
tive transport in two-dimensional flows^K^l For time- 
independent flows, these manifolds form separatrices that 
strictly separate distinct regions of the fluid (e.g. dif- 
ferent vortex cells.) For time-periodic flows, the sepa- 
ratrices split into distinct, but intersecting, stable and 
unstable manifolds, which form complicated heteroclinic 
(or homoclinic) "tangles"; advection between regions is 
governed through lobes formed by the intersecting mani- 
folds. More recently, invariant manifold techniques have 
been generalized to aperiodic and even turbulent flows 
using Lagrangian coherent structures and finite-time Ly- 
panunov exponent s^^. 

Here, we treat front propagation dynamics as a mod- 
ification of the above approach for passive advection. A 
point along the front is viewed as an active "particle" 
(or "swimmer") that evolves according to both fluid ad- 
vection and its own internal locomotion, fixed at a con- 
stant speed Vq (the "burning" speed) in the moving frame 
of the fluid. The corresponding dynamical system is 
thus three-dimensional, with two dimensions for the ap- 
position within the fluid and a third dimension for the 
orientation of the local front element. Initially, it may 
seem that the invariant manifolds defined for passive ad- 
vection are irrelevant for reactive front propagation, since 
the reaction can burn right through an advective sepa- 
ratrix. However — and this is the central point of this 
paper — the invariant manifolds of the advection dynam- 
ics can be suitably and naturally modified to the reaction 
front scenario; these are the burning invariant manifolds 
(BIMs). The BIMs reduce to traditional advective in- 
variant manifolds when vq — 0, but otherwise may differ 
dramatically from them. 

We give a brief overview of the physical role played by 




FIG. 2. (color online) BIMs in a time-periodic flow. An initial 
stimulation between the BIMs moves outside the BIMs during 
the forcing cycle (panel a) but returns by time T (panel b). 
A stimulation to the left of the BIM pair evolves through one 
BIM but not the other (panel c). Similarly for a stimulation 
on the right (panel d). 



BIMs. Figure [I] illustrates the BIMs (bold, red) for a 
time-independent flow, confined to a horizontal channel 
with counter-rotating vortices. In Fig.[]Ji, the reaction is 
catalyzed at the channel bottom between the BIMs, and 
then generates the sequence of reaction fronts, forming 
an upward propagating plume, ultimately rolled up by 
the vortices. The fronts are bounded by the BIMs on 
either side, with the arrows normal to the BIMs denot- 
ing the "burning" direction of each BIM. A BIM forms a 
barrier to fronts propagating in its burning direction (as 
in Fig. [IJi), but not to fronts impinging in the opposite 
direction, as illustrated in Figs. [l]o-d. In Figs. [TJ^-d the 
fronts are catalyzed at different stimulation points (small 
open circles) and approach the BIMs from different direc- 
tions. In each case, the fronts pass through a BIM when 
approaching opposite the BIM's burning direction but 
are stopped when approaching in the burning direction. 

In Fig. |2j the vortices oscillate periodically in the hor- 
izontal direction. The BIMs are shown at a single phase 
of the driving only. An initial stimulation between the 
BIMs in Fig. [2^ grows outward in a series of reaction 
fronts. During the first period, the front migrates out- 
side the BIM channel in Fig. [2ji, but returns to it exactly 
one driving period T after initialization (shown by the 
gray "burned" region in Fig.[2]3.) Figures]^ and[2]i show 
two additional stimulation points, and their subsequent 
evolution after two driving periods (in gray). In both 
cases, BIMs bound the front evolution in their burning 
direction, but not in the opposite direction. 

This paper is organized as follows. Section |II A| in- 
troduces the finite-dimensional dynamical system for a 
point along the front, and Sec. |IIB| discusses how ex- 
tended fronts are represented within this framework. 
Section |IID| compares this ODE approach to a grid- 
based computation. Front singularities and collisions 
are discussed in Sec. Ill El Section III Fl formalizes the 
observation that fronts cannot pass one another 
what we call the "no-passing" lemma. Section III de- 
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rives existence (Sec. Ill A) and stability (Sec. IIIC) re- 
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suits for fixed points of the front dynamics, focusing on 
time-independent flows. Linear flows are examined in 
Sec. |IIIB The existence of fixed points is constrained 



by a topological index theory developed in Sec. HID 



The heart of the paper is the introduction of BIMs in 
Sec. IV, which establishes that BIMs can (locally) be 
interpreted as fronts (Sec. IV A) that provide one-sided 
barriers to front propagation (Sec. IV B). Sections IV C 



and IV D discuss two mechanisms for circumventing these 



barriers. Section IV E describes how BIMs attract fronts, 
providing a mechanism to measure BIMs directly in the 
laboratory. Finally, Sec. [V| connects features of BIMs to 
the geometry of the simpler cases of either pure (passive) 
advection or pure (nonadvecting) front propagation. 



II. FRONT PROPAGATION IN AN ADVECTING 
MEDIUM: FUNDAMENTAL CONCEPTS 



We analyze the general problem of a front propagating 
through a two-dimensional medium flowing with velocity 
u(r, t), specified as a function of position r = (x,y) and 
time t. In the local co-moving frame of the medium, 
propagation is assumed to progress normal to the front 
at a speed vo (the burning speed), which is constant, 
independent of position, propagation direction, or the 
local front curvature^. Finally, the velocity field u(r, t) 
is assumed to be established independently from the front 
dynamics itself, i.e. the development of the front does not 
feedback to modify u. 

Our front propagation theory has been successfully ap- 
plied to the driven advection-reaction-diffusion (ARD) 
systems studied experimentally by Solomorl 18 * 2Q [ In these 
experiments, the reaction time scale is much faster than 
the advection time scale (i.e. large Damkohler number), 
meaning that the front is thin compared to the scale of 
the flow — the so-called "geometric optics" limit 28 ! In this 
limit, a fluid element is either fully reacted or unreacted; 
hence, the distribution of reaction products within the 
fluid is entirely determined by the boundary of the re- 
acted, or "burned" , region. Throughout this paper, such 
ARD systems form the prototype example, though our 
theoretical framework should lend insight into pattern 
formation in other systems that combine front propaga- 
tion and advection, including not only chemical reaction 
dynamics, but optics, phase-transitions, and even quan- 
tum phenomena. 



A. A three-dimensional ODE for the propagation of front 
elements in a two-dimensional flow 



A front element, i.e. a point on the boundary between 
burned and unburned regions, is described by both its 
position r = (x,y) and its unit normal n (Fig. [3|. Alter- 
natively, n can be replaced by the unit tangent g, defined 
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FIG. 3. Each reaction front element independently propa- 
gates forward under the ODEs j2| and J3}. 



here by a right-handed turn by tt/2 
g = Jn, where J 



-1 

1 



(1) 



Under the above assumptions, a front element evolves 

•,r;JHl 



r = u + «on 

= u - v lg, 
g = [g • (Vu)n]n 

= g • (Vu) - [g • (Vu)g] ; 



(2a) 
(2b) 
(2c) 
(2d) 



where dots denote time derivatives and Vu is the ve- 
locity gradient tensor with components (Vu)^ = 



i,j = x,y, where a comma denotes partial differenti- 
ation. The equivalence of Eqs. (2a) and (2b) and of 



Eqs. (2c) and (2d) is straightforward, the latter requir- 
ing |g]== 1. To justify Eqs. (2a) and (2d), first con- 
sider the nonburning case vq 



for which front elements 



should evolve under advection alone. Equation (2a) re- 
duces to the appropriate advective form r = u. Eq. (2d) 



remains unchanged; the first term generates the dynam- 
ics of advection on the (unnormalized) tangent vector, 
whereas the second enforces the normalization |g| = 1. 
[Equation (2d) readily yields d\g\ 2 /dt = 2g • g = when 
|g| = 1.] Considering the case vp ^ 0, the only neces- 



sary change to Eqs. (2a) and (2d) is to add the burning 
velocity vqii normal to the front in Eq. (2a). An appeal- 



ing physical analogy can be obtained by recognizing that 
Eq. ([2| also describes the motion of a thin rod-shaped 
swimmer aligned with g and propelling itself normal to 
its alignment direction at constant speed vo in the local 
fluid frame. 

Since |g| = 1 for all time, Eq. Q reduces to a three- 
dimensional ODE in (#,?/,#), where is the orientation 
angle between g and the positive x-axis (Fig. [3|, i.e. g = 
(cos#, sin0). 



y = u y 



■ vo sin#, 
vq cos#, 



= —v>x,x sine/ cost/ 

u VjX cos 2 + Uy^y 



■x,y 

sin cos 



sin 2 6>+ 



(3a) 
(3b) 

(3c) 



For incompressible fluids, u y , y = —u XyX , simplifying 
Eq. ( 3c ) . We emphasize that even though Eq. ^ is three- 
dimensional, the fluid flow itself is £u>o-dimensional. The 
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reduction from four [Eq. to three [Eq. ([3|] dimensions 
is useful for both numerical computations and many the- 
oretical arguments. 



B. Front compatibility criterion 

A front may either be viewed as a curve in the two- 
dimensional xy-sp&ce or in the three-dimensional xyO- 
space (or indeed even in the four-dimensional rg-space.) 
We clarify the relationship between these two represen- 
tations here. First, we parameterize fronts by their Eu- 
clidean length A, measured along the front in xy- space. 
By convention A increases in the +g direction; this is 
equivalent to keeping the burned region on your left while 
traversing the curve. Then any (differentiable) curve r(A) 
representing a front in xy-sp&ce lifts to a unique curve 
(r(A), #(A)) in x^-space, where 6(X) is the unique (mod- 
ulo 2tt) solution to 



(cos (9(A), sin (9(A)) = dr/dX. 



(4) 



On the other hand, a general curve (r(A),#(A)) in xyO- 
space is typically not the lift of a curve in xy-space — 
only those curves in xyO-space satisfying Eq. Q repre- 
sent fronts. We thus call Eq. Q the front compatibility 
criterion for a curve (r(A),#(A))T Geometrically, it is the 
requirement that at every point A of the curve, the two- 
dimensional vector dr/dX points in the direction g(A), 
defined by g(A) = (cos#(A), sin 0(A)). Thus the front 
compatibility criterion is equivalently written 



j dr dr 

dX dX 



(5) 



For curves in xy 6 -space, we reserve the term front for 
those curves satisfying the front compatibility criterion. 
To summarize, a front is represented either by a param- 
eterized curve in xy-sp^ce or a parameterized curve in 
xyO-sp&ce that satisfies the front compatibility criterion, 
and there is a one-to-one correspondence between these 
representations. 

If the front compatibility criterion [Eq. ([5|] holds along 
a curve (r(A, t), #(A, t)) at some initial time, then it must 
hold for all future times as the curve evolves. Physi- 
cally, this must be the case, since the time-evolution of a 
front must remain a front. To check that it follows from 
Eq. (§, assume first that Eq. ([5| holds at an initial time. 
Then, 



d 
dt 



dXj 
I ■ (Vu)njn • 
-g • (Vu)n -\ 



dr 

dX 



g.J 



(Vu) 



i dr 
dX 



v -i 



dX 



n • (Vu) T g = 0, 



(6) 



where the first equality follows from Eq. ([2| and the sec- 
ond from the fact that g • (dg/dX) = and Eqs. ([I]) and 
Q. Thus, if Eq. ([5| holds at an initial time, it holds for 
all future times. 
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FIG. 4. a) An initial circular burned region at t — 0.25T, 
with the black front surrounding it; Panels b-d show the for- 
ward evolution at times t = 0.69T, 0.97T, 1.16T. The black 
front is evolved using Eq. |3| with vq — 0.2 and u given by 
Eq. Q, with b = 0.25, and uo — 2. The gray burned domain 
is evolved according to a grid-based, split-step algorithm, us- 
ing 1668x1668 grid points within the x2/-domain shown (full 
integration domain is larger), with time step At = 27t/(32cj). 



C. Example Fluid Flow 

Equations Q and ([3| are valid for an arbitrary flow 
u(r , t). For illustrative purposes (see Figs. 0§ [4) [10) and 
pT]), we use the following incompressible velocity field 



u x {x, y,t) = + sin(7r[x + bsm(ujt)]) cos(7ry), 
u y (x, y,i) = — cos(7r[x + bsm(uji)]) smf^ny). 



(7) 



This flow represents an array of alternating vortices, pe- 
riodically driven with frequency uo and amplitude b. It 
has been exten sively used to su ccessfully model driven 
laboratory flowP 



D. Two approaches to front propagation: curve-based 
versus grid-based 



Under the assumptions laid out in Sec. |II A[ an entire 
front can be evolved by independently evolving individ- 
ual front elements. Numerically, this requires a sufficient 
density of points along the front to resolve it to a desired 
accuracy. As the front evolves, the separation between 
some of the points typically increases, requiring the ad- 
dition of new interpolated points to maintain sufficient 
density. We use an algorithm similar to Ref. [40j 
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Refs.l38 l l4T | l42l used a Euclidean-grid-based approach 
to evolve the entire burned region itself. In such an ap- 
proach, each grid point is recorded as either burned or 
unburned at a given time. A split-step scheme, consist- 
ing of a pure burning step and a pure advection step, 
is used to update each gridpoint in time. (For details, 
see Refs.[38j[4TJ|42j) In contrast, our approach focuses 
solely on the front itself. 

Fig. [4] compares the above two algorithms for front evo- 
lution. Fig. [4^i shows an initial burned circular region. 
The black curve denotes the front (x(A), 2/(A), 0(A)) that 
forms the initial condition (at t = 0.25T, with T = 2tt/uj) 
for evolution via Eqs. (|3| and ([7]). A corresponding grid- 
based computation is also performed, with initial condi- 
tion consisting of the burned grid points inside the circle 
(shaded gray.) Fig. shows the independent evolution 
of both the black circle [via Eq. p|] and the gray grid 
points [via the split-step algorithmfup tot = 0.69T. The 
two techniques agree well, with the black curve marking 
the boundary of the gray region, thereby providing added 
numerical confirmation for the validity of Eq. The 
tiny visible discrepancy is due to the finite grid resolu- 
tion. 



E. Front intersections and cusps 

Fig. [4)3 continues the evolution to t = 0.97T. Part 
of the evolved front now lies within the burned (gray) 
region. While we refer to the entire evolved curve (ei- 
ther in xy-sp&ce or x?/0-space) as the front, we refer to 
those segments of the front separating the burned from 
unburned regions as the bounding front; those intervals of 
the front not in the bounding front thus lie interior to the 
burned region. Note that the bounding front need not be 
connected, but can have disconnected segments defining 
voids, as seen at time t = 1.16T (Fig.[4]i). Figures^ and 
[4]i continue the excellent agreement between the evolu- 
tion algorithms. Note that even though the front may 
intersect itself when projected into the xy-plane, it can- 
not self intersect in the full x?/0-space, as mandated by 
the uniqueness theorem for solutions of ODEs. 

Figures ^ and [4]i illustrate two mechanisms by which 
a front may penetrate the interior of the burned region. 
In the transition from Fig. |4Jd to Fig. [4)3, a so-called 
swallowt ail ca tastrophe is formed by the ^-projection of 
the frontP^ES the characteristic swallowtail shape, with 
its two cusps, is seen interior to the burned region in 
Fig. [4)3. Despite singularities in its ^-projection, the 
front in xyO- space remains smooth throughout. In the 
transition from Fig.[4]3 to Fig.|4]i, fronts penetrate the in- 
terior through the head-on collision of distinct locations 
along the front. Both mechanisms for producing front 
intersections in t he x ?/-plane are characteristic of front 
evolution in opticP^^, for example, and do not require 
the presence of advection. 



F. Front no-passing lemma 

Consider two fronts, with one front trailing the other. 
Since both fronts have the same speed in the local fluid 
frame, it is clear that the trailing front can not catch up 
to and pass the leading front. This no-passing result will 
be central to understanding why BIMs form one-sided 
barriers to front evolution in Sec. IIVI 

The intuitive no-passing statement above can be recast 
more precisely as follows. Let Co and Cq be the leading 
and trailing fronts, and assume that they are represented 
by (piecewise smooth) nonintersecting closed loops in xy- 
space, with Cq nested inside Co (Fig. |5^i). Without loss 
of generality, we assume that both are burning outward, 
meaning that their interior domains are burned. (The 
case of both burning inward, with their interior domains 
unburned, can be handled analogously.) The fronts may 
extend to infinity or coincide with a boundary wall for 
part of their length, so long as the burned region of Cq 
is entirely included within the burned region of Co. As 
Co and Cq evolve forward for some time t, according to 
Eq. (|3|, they produce new fronts C t and C[. In general, 
Ct and C[ might self-intersect and intersect one another 
(in the x7/-plane), as shown in Fig. [5J1. However, de- 
spite these intersections, the bounding fronts Ct and C[ 
(in bold) will not intersect and will retain their original 
nested ordering, with C' t lying within the burned region 
defined by Ct- We call this the no-passing lemma. It 
follows from the simple physical observation that the in- 
nermost burned region must always remain a subset of 
the outermost burned region. 

The no-passing lemma can be recast as a local formula- 
tion in terms of non-closed segments of the fronts C t and 
C' t . Suppose (rj, 0' t ) is a front element (of curve C' t ) whose 
position y[ coincides with a point r t on the front Ct, but 
whose 0' t does not agree with the local orientation of Ct. 
Then a short time e earlier, the front element (r^_ e , 0' t _ e ) 
must lie in front of C t - e (Figs. [5]o and [5^). Said another 
way, the front element (rj,0£) can not catch up to and 
pass Ct from behind. Of course, there is nothing pre- 
venting (r^,0Q from colliding with C t from the front, as 
in Figs. [5)3 and [5)3. Formally, the local no-passing formu- 
lation may be recast as the inequality 

r' • n < r • n, (8) 

where r' and r are the velocities of the two front elements 
at the time their ^-positions coincide and n is the unit 
normal to C^. Equation ([8| follows directly from r = 
u(r) + von [Eq. ([2^l)], r' = u(r) + von' , and n • n' < 1. 
Finally, note that the roles of (r£,0£) and (v tl t ) can be 
reversed, so that (r t , t ) can not catch up to and pass C' t 
either; mathematically, r • n' < r' • n'. 

III. BURNING FIXED POINTS 

By a burning fixed point, we mean either an equilibrium 
point of the differential equation ^ [or ([3|], for time- 
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FIG. 5. a) The fronts Co and Cq evolve into Ct and C' t . The 
bounding fronts C t and C' t (in bold) have the same nested 
ordering as Co and Cq. Panels b and c illustrate a front 
element (rj,0') (of C' tl not shown) approaching Ct from the 
front and colliding with it. 



independent flows, or a fixed point of the corresponding 
Poincare map, for time-periodic flows. For the remain- 
der of Sec. |III[ we restrict our focus to time-independent 
flows, for which we give a precise characterization of the 
existence and stability of burning fixed points. 



A. Existence criteria for burning fixed points in 
time-independent flows 

At a burning fixed point in a time-independent flow, 
Eq. (2a) implies vqu = — u, so that the fluid velocity is 



exactly balanced by the burning velocity. Eq. ( 2c ) on the 
other hand implies = n • g = g • (Vu)n. 

Theorem 1 A burning fixed point (r, g) occurs if and 
only if the following two conditions are met. 



(0 

(ii) 



v n : 
g • (Vu)n = 



-u, 



0. 



Condition (ii) is equivalent to 

(ii) (Vu) T g = /ig, 
for some (necessarily real) eigenvalue \i. 



(9) 
(10) 



(11) 



A n 



g 



C 



u 




A*>0 



fi<0 



FIG. 6. Illustration of flows with \i > and \i < about a 
burning fixed point, at which u is normal to the set C defined 
by |u| = vq. The burning direction n points opposite u. 



Together, Eqs. (|9 
is tangent to the level set 



and (10) imply g • V|u| 2 = 0, i.e. g 
= vq (assuming V|u| 2 =^ 0). 



Alternatively, n is normal to the level set, so that by 
Eq. ([9| u is also normal to the level set. 

Theorem 2 Define the level set C in xy -space by |u| = 
vq. Assuming V|u| 2 7^ on C, the set of burning fixed 
points consists of exactly those points (r, g) at which r is 
on C and u(r) points normal to C. The burning direction 
n = — Jg points opposite u. 

Figure [6] illustrates Theorem [2] for /i > and /i < 0. 



B. Burning fixed points for linear fluid flows 

The case of linear flows lends insight into the existence 
and local structure of burning fixed points. 



1. Hyperbolic flows 



The hyperbolic flow 



-Ax, 



(12) 



with strength A > 0, has a unique hyperbolic advective 
fixed point at the origin. This flow generates the burning 
dynamics 



x = — Ax + sin#, 
y = +Ay — vq cos 
6 = +2A cos <9 sin <9. 



(13a) 
(13b) 
(13c) 



Thus, any fixed point (x* : y* : 6*) must have 6* = 
0,7r,7r/2,37r/2. Each 9* then yields a unique value of 
x* and summarized in Table [IJ The advective hyper- 
bolic fixed point thus produces four burning fixed points 
as deviates from 0. 

Figure [7^ shows the four burning fixed points. Consis- 
tent with Theorem [2j these four points lie on the circle 
A\r\ = |u| = i>o, bounding the shaded region |u| < vq. 
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BFP (x*,y*,6g 


Eigenvalues 


Stability 


ID BIM dir 


1 


(0, VAO) 




SUU 


X 


2 


(0,-vo/A,tt) 


-A,+A,+2A 


SUU 


X 


3 


(VA0,7r/2) 


-A,+A,-2A 


ssu 


y 


4 


(-wo/A,0,37r/2) 


-A,+A,-2A 


ssu 


y 



TABLE I. Burning fixed points for hyperbolic flow Eq. ( 12 ) 




FIG. 7. (color online) a) The four burning fixed points (filled 
dots) in a hyperbolic flow lie on the circle C (the level set |u| = 
vo) which surrounds the advective fixed point (open dot). The 
shaded region satisfies |u| < vo. The burning velocities (large 
arrows) of these four points are normal to the circle, pointing 
opposite the flow. The left and right points have stability 
SSU, and the vertical lines passing through them are their 
unstable BIMs. The triangular (unbarbed) arrows denote the 
unstable direction along the BIMs. b) An elliptic flow has no 
burning fixed points since the flow lines are never normal to 
the ellipse C bounding the shaded region. 



Furthermore, these are exactly the four points where the 
hyperbolic flow is normal to the circle. 



The linearization of Eq. (13) yields 



d(x,y,0) 



-A vo cos 
A vq sin 
2 A cos 20 



(14) 



The eigenvalues of L at the burning fixed points 
are (—A, A, ±2 A), summarized in Table |lj Each fixed 
point is hyperbolic with either stability stable-stable- 
unstable (SSU) or stable-unstable-unstable (SUU). Each 
fixed point with SSU (SUU) stability thus has a one- 
dimensional unstable (stable) invariant manifold, which 
is a burning invariant manifold (BIM). These manifolds 
have eigendirections (y, 0) (for SSU) or (x, 0) (for SUU). 
In both cases, the ^-projection of the eigendirection is 
aligned with the orientation g of the fixed point. 

Focusing on the two burning fixed points on the 
x-axis, their unstable BIMs are the straight lines 
{(±vq/A, y, 7r/2) I — oo < y < oo}, which project to the 
two vertical lines in Fig. [7^i. Note that these BIMs sat- 
isfy the front compatibility condition Eq. Q, a point we 
return to in Sec. IIV A[ 



2. Elliptic flows 

The elliptic flow 

u x y/A, u y -Ax, 



(15) 



with A > 0, has a unique elliptic advective fixed point at 
the origin. It generates the burning dynamics, 

x = y/A + vo sin (16a) 
y = — Ax — vq cos#, (16b) 
= -(sin 2 (9) /A- A cos 2 (9, (16c) 

from which we see that = has no solutions, and hence 
there are no burning fixed points in a purely elliptic flow. 
This can also be seen as a consequence of Theorem [2j as 
illustrated in Fig. [7}}. Any burning fixed point must lie 
on the elliptic level set [(Ax) 2 + (y/A) 2 ] 1 / 2 = |u| = v 
bounding the shaded region. However, the flow is never 
normal to this set, and hence there are no burning fixed 
points. 



C. Stability of burning fixed points for general 
time-independent flows 

For the main results on the stability of burning fixed 
points, the reader may skip ahead to Theorems [3] and [4] 
The linearization of Eq. ([2| yields the following four-by- 
four matrix, decomposed into two- by- two blocks. 



X 



where 



dr/dr 


dr/dg 


dg/dr 


dg/dg 



' (Vuf 


-v J 


n (8) a 


D 



(17) 



D = (g (g> n + n <g> g) [g • (Vu)n] ■ 

n ® n [Tr(Vu) - 2g • (Vu)g] , 
a = g • VV(n • u). 



(18) 
(19) 



Here, a tensor product b(g)c has components (b <g> c)ij = 
biCj and a$ = gknjUj,kn wnere repeated indices are 
summed according to the Einstein convention. 

We now specialize to linearization about a burn- 
ing fixed point for time-independent flows, and assume 
V|u| 2 ^ 0. Equations (T 



D 



n ( 



lepen 
9j) - ([n) imply 

> n [Tr(Vu) - 2/j] . 



(20) 



Now, since |g| = 1 is preserved by the dynamics, we 
need only consider variations of g in the direction n; i.e. 
we consider the three-dimensional space of vectors of the 
form (v,cn), with v £ R 2 and c £ R, which is invariant 
under X. Restricting X to this space yields the three-by- 
three matrix 



Y 




-^og 
Tr(Vu) - 2/i 



(21) 
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Equation (11) implies that the two-dimensional space of 
vectors of the form (ag, c), with a, c G M, is invariant 
under Y. Restricting Y to this space yields the two-by- 
two matrix 



li -v 
a • g Tr(Vu) - 2/i 



(22) 



The eigenvalue of Y whose eigenvector does not have the 
form (ag, c) is thus 

A = Tr(Y) - Tr(Z) = Tr(Vu) - /i = v, (23) 

where v is the remaining (necessarily real) eigenvalue of 
(Vu) T . The remaining two eigenvalues of Y are thus the 
two eigenvalues of Z given by 



Ah 



\{y ± V" 2 -4[v a.g + /i(i/-/i)]) . 



(24) 



This expression can be put in a simpler form if we first 
consider the components of u in the local frame of C (the 
level set |u| = i?o), defined by e[ and e^, the unit normal 
and tangent vectors to C, respectively. The components 
of u in the (e^e^) frame are denoted u' = (u^u^). We 
parameterize C by the Euclidean distance A measured 
along C. At a burning fixed point, the derivative of u 
along C is thus du/dA, which by Eq. ( [TT] ) equals /xg, 
taking A to increase in the direction g. The quantity 
fi/vo = fi/\u\ is thus a kind of rotation rate of u along 
C. Similarly, we may differentiate u in the (e' l5 e^) frame, 
i.e. du'/dX. Since \w'\ = vo is constant along C, du' /dX 
must point orthogonal to u' . At a burning fixed point, g 
is orthogonal to u, and hence 



du'/dA = //g' 



(25) 



where g' = (0, 1) is the components of g in the (e'^e^) 
frame. The (real) quantity [i'/vq is thus also the rota- 
tion rate of u along C, but viewed in the (e^e^) frame. 
Appendix [A] shows 

[/ = fj, + VqK, (26) 

where n is the (signed) curvature of C, which equals (see 
Appendix |A| ) 



v a • g - /T 



(27) 



at a burning fixed point. Combining Eqs. (26) and (27) 
yields 



= - boa • g + [i[y - , 



(28) 



which can be used to transform Eq. (24), as below. 



Theorem 3 For a time-independent flow u, the eigen- 
values about a burning fixed point are 



An = v. 



4z//i' 



(29) 
(30) 



where v is the eigenvalue o/(Vu) T not given by Eq. 
and // is given by Eq. 



\2by . The linear stability of a 
burning fixed point is thus determined by the signs of v 
and ji' according to the following table. 





v < v > 




> 


ssu uuu 


// 


< 


sss suu 



In the above, v determines the stability of two eigen- 
values and // the stability of one. In either the SSU or 
SUU case, all eigenvalues are real. However, in the UUU 
or SSS case, two eigenvalues are complex if A/a'/v > 1. 
In Sect. [IV] we shall pay particular attention to the SSU 
case. For incompressible flows v = —fi, and Eq. (28) 
reduces to 



(i 



von ■ gl 



(31) 



and Theorem [3] specializes to the following. 

Theorem 4 For a time-independent, incompressible 
flow u, the eigenvalues about a burning fixed point are 



A 
A± 



-/i, 



4/i/V 



(32) 
(33) 



11) and (25). The 



where \i and \j! are given by Eqs. 
linear stability of a burning fixed point is thus determined 
by the signs of \i and ji' according to the following table. 



Ii > 11 < 



li' >0 
li' < 



SSU uuu 
sss suu 



For a linear flow (a = 0), Eqs. (32) and (33) reduce 
to the eigenvalues in Table [IJ More generally, for either 
sufficiently small or sufficiently linear flow (a small), 
Eq. (31) shows that // has the same sign as /i, and thus 



a burning fixed point has one of only two stability types: 
SSU (ii > 0) or SUU (ii < 0). Referring to Fig. § defo- 
cusing flows yield SSU and focusing flows SUU, as verified 
by the fixed points in Fig. [7k. 



D. Topological index theory for the existence of burning 
fixed points 

The number and stability type of burning fixed points 
are constrained by a topological index theory similar 
to the Poincare index theory for fixed points in two- 
dimensional flowP'. Following Theorem |2j we consider 
a closed loop C upon which |u| = ^o- Since |u| ^ on 
this loop, the winding of the orientation of u can be used 
to define two topological indices. The first is the usual 
Poincare index^. 

np(C) = the number of counterclockwise (ccw) rota- 
tions made by u under one ccw circuit of C. 

(34) 








FIG. 8. Four flows illustrating Poincare and burning indices. 



FIG. 9. The burning index ±1 is given for each velocity con- 
figuration. The shaded region satisfies |u| < vq. 



For each rotation made by u', u' is perpendicular to 
C twice, once on either side of C. Thus, there are at 
least \ub\ outward-burning fixed points on C and \np\ 
inward-burning fixed points. This result can be refined 
by defining a burning index for each burning fixed point 



Fig. [8] illustrates rip for four example flows. The reader 
is encouraged to confirm the value of rip cited in each 
case. Following standard practice, a Poincare index is 
also defined for each advective fixed point z as follows. 



n P (z) 



+1 z elliptic 
— 1 z hyperbolic 



(35) 



The Poincare index of the loop C is then equal to the sum 
of the Poincare indices of all the enclosed fixed points 

np(C)= n P ( Z ) = N e - N h , (36) 

fixed points z 

where N e and Nh are the number of elliptic and hyper- 
bolic advective fixed points, respectively. The reader can 
confirm this relationship in Fig. [8] 

Alternatively, by considering the components u' in the 
local frame (e^e^) of C, we define the burning index as 
follows. 

tib(C) = the number of ccw rotations made by u' under 
one ccw circuit of C. 

(37) 

The two indices are related by 



n B (C)=n p (C) 



(38) 



since the vector u in the original laboratory frame un- 
dergoes one additional rotation relative to the vector u' 
in the local frame of C. 



, . I +1 u rotates ccw at z, traversing C ccw 
n B {z) < 

1—1 u rotates cw at z, traversing C ccw 

(39) 

This definition is illustrated by the four cases in Fig. [9] 
These four cases also make it clear that the burning in- 
dex of a fixed point can be re-expressed in terms of \jl 
[Eq. ([25b] as 



n B (z) 




and u(z) points inward to C 
(outward-burning) 



and u(z) points outward to C 
(inward-burning) 

(40) 

Analogous to the Poincare index, the burning index 
npiC) of the curve C is the sum of the burning indices 
of either all the outward-burning fixed points or all the 
inward-burning fixed points located on C. 



n B (C) 



E 

: BFI 

E 



n B (z) = N° ut - N% 



out 



out BFPs z 



n B (z) = N^-Ni! 



(41) 
(42) 



in BFPs z 



where 7V° ut and N° ut are the number of outward-burning 
fixed points with \j! > and // < 0, respectively, and 
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and are the number of inward-burning fixed points 
with // > and // < 0, respectively. Finally, the number 
and type of burning fixed points on C can be directly 
related to the advective fixed points interior to C through 
Eqs. d36b and d38k. 



For an incompressible flow, \i must be positive when 
u points inward and negative when u points outward. 
Using the table in Theorem [4j Eq. (40) then becomes 



n B (z) = < 



+1 stability SSS and outward-burning 

— 1 stability SSU and outward-burning 

+1 stability UUU and inward-burning 

— 1 stability SUU and inward-burning 



Similarly, Eq. (43) becomes 



(44) 



Nsss ~ N °ssu = Nuuu - Nlhu = N e -N h -l, (45) 

where Nxxx measures the number of inward/outward- 
burning fixed points of the given stability type. This 
equation is confirmed in Fig. |§J where all burning fixed 
points either have stability SSU or SUU and hence the 
number of outward- or inward-burning fixed points sim- 
ply equals 1 — N e + Nh = — n#(C). 



IV. BURNING INVARIANT MANIFOLDS (BIMS): 
ONE-SIDED BARRIERS TO FRONT PROPAGATION 

Having examined the existence and stability of burning 
fixed points in Sect. |III[ we now focus on the invariant 
manifolds attached to them. These manifolds are invari- 
ant either under the differential equation ([2| [or Q], in 
the time-independent case, or under the Poincare map, in 
the time-periodic case; a stable/unstable invariant man- 
ifold consists of all points that converge upon a burning 
fixed point in forward/backward time™. As noted above, 
we call these stable and unstable manifolds burning in- 
variant manifolds (BIMs) to distinguish them from the 
traditional invariant manifolds for advection. Due to the 
special role played by one-dimensional invariant mani- 
folds, all BIMs subsequently discussed are assumed to be 
one-dimensional, unstable manifolds, attached to a fixed 
point of stability SSU, as shown in Figs. [I] and [2] for time- 
independent and time-periodic flows. (One-dimensional, 
stable BIMs are similarly analyzed by running time back- 
wards; they will be relevant in a future publication for 
developing a theory of lobe dynamics for front propaga- 
tion.) The central point of our theoretical development 
is that these BIMs form one-sided barriers to front prop- 
agation (Figs, [I] and [2J. 



A. BIMs are fronts 

Though most one-dimensional curves in the three- 
dimensional x^-space are not fronts, an unstable (one- 
dimensional) BIM attached to an SSU fixed point does 
satisfy the front compatibility criterion, i.e. Eq. Q or 
([5]), as established here. Consider first a short line seg- 
ment y*, #*) + (Acos(0), Asin(#), 0), parameterized 
by — e < A < e and passing through a burning fixed 
point (#*,?/*,#*) of stability SSU. This line clearly satis- 
fies Eq. Q. Furthermore, due to the two stable dimen- 
sions of the fixed point, as this curve evolves in time, 
it converges upon the BIM, approximating progressively 
longer segments of the BIM ™ . Since the front com- 
patibility criterion is preserved over time, the evolving 
curve remains a front, and hence the BIM to which it 
converges satisfies the front compatibility criterion every- 
where along its length. This argument applies for both 
time-independent and time-periodic flows. 



B. BIMs are one-sided barriers 

As pointed out in Figs, [I] and |2j BIMs are one-sided 
barriers, preventing fronts from passing in one direction, 
but not the other. We can now explain this fact as fol- 
lows. Since B IMs are fronts, they satisfy the no-passing 
lemma (Sect. II F). The local formulation of this lemma 



then implies that no impinging front can burn past any 
point of the BIM in the same direction that the BIM 
itself is burning. 

In some flows, a BIM can partition the fluid into two 
distinct regions, one entirely burned and one entirely un- 
burned. The BIM then prevents the burned region from 
penetrating the unburned region. This is seen, for exam- 
ple, in the hyperbolic flow in Fig. [7^i, in which the BIMs 
form two vertical lines. If the entire vertical strip be- 
tween the two BIMs were burned, the region outside the 
strip would remain unburned forever; consequently, a re- 
action catalyzed anywhere between the BIMs will never 
burn outside the central strip. Alternatively, the BIMs 
in Figs, [l] and [2] do not confine the burned material to 
a bounded region forever; in these examples, fronts will 
eventually reach the other side of any given BIM. This 
can happen in two ways, discussed below. 



C. Cusps form openings in the barrier 



The 
in Figs 
on the 
burning 
BIM is 
See Fig 
toward 
front, 
toward 



^-projection of a BIM may form cusps, as seen 
.[I] and [2] A cusp in a BIM has a profound impact 
bounding nature of the BIM. This is because the 
g direction (and hence bounding direction) of a 
reversed at a cusp, relative to an impinging front. 



10 i. To the right of the cusp, the BIM is burning 



the upper left, acting as a barrier to the impinging 
However, left of the cusp, the BIM is burning 
the lower right. This flip in the burning direction 
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front 



FIG. 10. (color online) a) A BIM cusp, b) Zoom out of cusp 
with fronts wrapping around, c) Fronts spiraling and sliding 
past the cusp, d) A BIM with swallowtail structure. 



results from the continuity of 6 along the curve in the 



full xyO-spdice — note that the arrows in Fig. 10 1 rotate 
smoothly along the BIM. Thus, left of the cusp, the BIM 
does not act as a barrier to the front below; the front will 
simply pass through the left segment, as in Fig.lTofc. Left 
of the cusp, the front spirals clockwise around the cusp 
and forms a swallowtail singularity. The front then slides 
past the cusp itself, moving toward the lower right and 
crossing over the right BIM segment from above. The 
critical observation is that the cusp provides an opening, 
through which the front can pass around the bounding 
segment of the BIM, and then collide with it. For time- 
independent flows, this is the only way for a front to reach 
the other side of a BIM. 

Figure [To} ) shows a zoomed out view of the BIM and 
an initial front impinging from the left. Note that the 
evolving front passes right through the continuation of 
the BIM segment shown left of the cusp in Fig. 10 1, ow- 



ing to the orientation of its burning direction. However, 
although a cusp creates an opening in the bounding na- 
ture of a BIM, the entire remainder of the BIM is not 



necessarily irrelevant. Figure 10 i shows an example of 
how a second cusp again reverses the burning direction^. 
Though the front may burn around each of the two cusps, 
the (bold) intervals still form a barrier to further front 
propagation. 



D. Moving past BIMs in time-periodic flows 

For a time-periodic flow, a front may reach the other 
side of a BIM by a second mechanism illustrated in 
Figure [TT^ i shows an initial segment of a left- 
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Fig. 11 



burning BIM (bold, red) that transects the fluid channel. 



FIG. 11. (color online) The mechanism by which a burned 
region (gray) may move past an initial BIM segment (red) 
under time-periodic driving. 



A gray burned region is immediately right of the BIM. As 
the system evolves over one driving period, the BIM seg- 
ment stretches and folds, ultimately crossing over itself 
(in the ^-projection). By the final frame, both the BIM 
and the front extend left of the original BIM segment 
(bold). While this may appear to violate the no-passing 
lemma, at no point does the burned region cross over any 
point of the BIM in the burning direction of the BIM. 
Furthermore, at each time step, the leftmost transecting 
segment of the BIM bounds the burned region. By the 
final frame, the burned region has moved past the origi- 
nal BIM segment because of the stretching and folding of 
the BIM itself. This mechanism is analogous to that of 
lobe dynamics (or turnstiles) for the advec tive tran sport 
of passive impurities in time-periodic flow d 9 * 11 * 12 * 48 *. 



E. Convergence to BIMs and experimental measurement 

Fronts are not only bounded by BIMs but converge 
upon them. This is due to the transverse stability of a 
BIM in the neighborho od of i ts SSU fixed point (a fact 
already utilized in Sec. IV A). If a front has one front 
element that converges on a burning fixed pointP^, nearby 
front elements will evolve arbitrarily close to any given 
point of the BIM over time. This convergence is clearly 
witnessed in Figs. [T] and [2Jd, d. 

Front convergence suggests a protocol to measure 
BIMs in the laboratory. For the time-independent case, 
first catalyze a reaction front sufficiently near the burn- 
ing fixed point. (The advective fixed point is often a good 
choice.) Then observe the limiting behavior of the evolv- 
ing front, which identifies the BIM location, as simulated 
in Fig. Points on the BIM nearest the fixed point 
will approach their limiting value earlier than points fur- 
ther along the BIM, and hence can be measured with 
greater accurac}^. This protocol has been successfully 
implemented in the Solomon lab, as discussed in the com- 
panion article 20 and Ref. US 

For the time-periodic case, a similar protocol exists, 
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Front Propagation Effects 

• Growth in the total area of the impurity region 

• Formation of cusps, swallowtails, and self-intersections at 
the impurity boundary (and within a BIM) 

• Directionality (i.e. existence of a burning direction) for the 
impurity boundary (and for a BIM) 

Advection Effects 

• Existence of fixed points and invariant manifolds (BIMs) 

• Bounding property of invariant manifolds and use in trans- 
port 

• Convergence of impurity boundary to invariant manifolds 

• Lobe dynamics approach to transport of impurity region 

TABLE II. Comparison of front versus advection effects 



passive tracers or an active chemical medium) within a 
"background" medium. The front dynamics, and the 
BIMs in particular, reflect properties unique to one limit 
or the other, as summarized in Table [HI For example, the 
formation of cusps in a BIM reflects behavior typical in 
front propagation, but does not occur for invariant man- 
ifolds used in advection. Similarly, BIMs have a burning 
direction, which is natural for fronts, but absent for ad- 
vective manifolds. On the other hand, the very existence 
of BIMs, their use as bounding curves, and their appli- 
cation to transport reflect similarities with the advection 
problem but are absent from studies of front propagation 
in homogeneous stationary media. 
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with one critical change. Care must be taken to record 
the converging front at a constant phase of the driv- 
ing (say t = 0). Hence, multiple experimental runs are 
needed in which the reaction is stimulated at incremen- 
tally earlier times t < 0, and then recorded at the com- 
mon time t = 0. See Ref. 18 for further details and an 
experimental implementation. 

BIMs can also be computed directly from the dynamics 
Eq. (|3| if the flow u(r, t) is known. However, the abil- 
ity to measure BIMs directly in experimental flows frees 
one from the requirement of either theoretically model- 
ing u(r, t) or experimentally measuring u(r,t), e.g. from 
particle-tracking data. 

As a front is converging upon a given point r of a BIM, 
another interval of the front may have taken a different 
route to collide with the BIM at r from the other side 



(as discussed in Sects. IV C and IV D). Once this hap- 
pens, one can no longer image the convergence upon r, 
since the entire neighborhood of r is burned. Thus, the 
accuracy by which the BIM can be measured is limited 
by the time it takes another interval of the front to col- 
lide with the BIM. The best strategy then is to initialize 
the front as close to the burning fixed point as possible, 
allowing the most time to approach any given point on 
the BIM before the convergence process is wiped out. An 
enhanced strategy, permitting a more detailed measure- 
ment of the BIMs, would be to also suppress, or "erase" , 
the second front before it collides with the point r be- 
ing measured. This could be done optically for certain 
photosensitive reactions. 



V. CONCLUDING REMARKS 

A useful framework for organizing the properties of 
front propagation in a flowing medium [Eq. ([2|] is to 
consider the two natural limiting cases: u = (front 
propagation in a stationary medium) and = (pas- 
sive advection). In both limits, one is concerned with 
the change in geometry of an "impurity" region (either 
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Appendix A: Derivations 
1. Derivation of Eq. ( [26] ) 

The (signed) curvature of a planar curve C equals 



dA ' 



(Al) 



where and e' 2 are the norma l and tangent unit vectors 
to the curve, as in Sec. Ill C| (Here, k < implies 
points toward the center of curvature.) We denote R(A) 
as the rotation matrix connecting the (e^ej?) frame to 
the laboratory frame (ei,e2) = (x, y), i.e. 

eJ(A) = R(A)e z , z = 1,2. (A2) 

The components of a vector then transform using R T , e.g. 

u'(A) = R T (A)u(A). (A3) 

Specializing to a burning fixed point, at which n = 



and g = e' 2 (Theorem |2j, Eqs. |AT} and flA2| yield 



dR. 

dX 



R T e[=g 



dR 
d\' 



R 1 n. 



(A4) 



From du/dX = /ig, we then have 



M = g 



du 

dX 



dR^ T ^du' 
dX dX 



dR^ T 

dX 



VR g ' 



-vqk + h', (A5) 



where the second equality follows from Eq. (A3), the 
third from Eqs. Q and (25), and the last from Eqs. (A4) 
and the fact g = Rg / . This yields Eq. (26). 
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2. Derivation of Eq. ([27) 

The (signed) curvature of the level set /(r) 
general function / is 

_ e' 2 ■ (VV/)e' 2 
K ~ V/-ei • 



for a 



(A6) 



Taking / = |u| 2 — defines the level set C. Since e[ — n 
and e' 2 = g at a burning fixed point, the denominator of 
Eq. (|A6|) becomes 



(V|u| 2 ) - n = -2v n- (Vu)n 

= -2^(Tr(Vu)-g-(Vu)g) 



(A7) 



where the first equality follows from Eq. ^ and the final 
equality from Eq. (11) and the fact that Tr(Vu) = /i + zA 
Similarly, a straightforward applic atio n of Eqs. (Jo]) and 
(11) reveals the numerator of Eq. (|A6[) to be 



g-(VV|u| 2 )g = 2/i 2 -2^a-g. 



(A8) 



Eq. (27) then follows from Eqs. (A6) - (A8). 
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